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ABSTRACT 

We present results of the multi wavelength campaign on the TeV blazar Mkn 501 
performed in 2006 July, including MAGIC for the very-high-energy (VHE) 7-ray band 
and Suzaku for the X-ray band. A VHE 7-ray signal was clearly detected with an 
average flux above 200 GeV of ~20% of the Crab Nebula flux, which indicates a low 
state of source activity in this energy range. No significant variability has been found 
during the campaign. The VHE 7-ray spectrum can be described by a simple power- 
law from 80 GeV to 2 TeV with a photon index of 2.8 ±0.1, which corresponds to 
one of the steepest photon indices observed in this energy range so far for this object. 
The X-ray spectrum covers a wide range from 0.6 to 40 keV, and is well described by 
a broken power law, with photon indices of 2.257 ± 0.004 and 2.420 ± 0.012 below 
and above the break energy of 3.24+q }2 keV. No apparent high-energy cut off is seen 
above the break energy. Although an increase of the flux of about 50 % is observed 
in the X-ray band within the observation, the data indicate a consistently low state of 
activity for this source. Time-resolved spectra show an evidence for spectral hardening 
with a flux level. A homogeneous one-zone synchrotron self-Compton (SSC) model 
can adequately describe the spectral energy distribution (SED) from the X-ray to the 
VHE 7-ray bands with a magnetic field intensity B = 0.313 G and a Doppler beaming 
factor 5 = 20, which are similar to the values in the past multiwavelength campaigns 
in high states. Based on our SSC parameters derived for the low state, we are able 
to reproduce the SED of the high state by just changing the Lorentz factor of the 
electrons corresponding to the break energy in the primary electron spectrum. This 
suggests that the variation of the injected electron population in the jet is responsible 
for the observed low-high state variation of the SED. 

Subject headings: BL Lacertae objects: individual (Markarian 501) - galaxies: jets - 
gamma rays: observations - X-rays: galaxies 



1. Introduction 

Blazars are radio-loud active galactic nuclei (AGNs) viewed at small angles between the jet 
axis and our line of sight. The bulk relativistic motion of the emitting plasma causes the radiation 

"•"deceased 
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physics and Cosmology, SLAC National Accelerator Laboratory, CA, 94025, USA 
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to be beamed in a forward direction, making the variability appear more rapid and t he luminos- 
ity appear higher tha n in the rest frame due to the relativistic beaming effect (e.g., iReesI Il966 
Ghiselliniet all 1993b . 



Blazars with only weak or entirely absent emission lines in the optical band are classified 



as BL Lacertae Objects (BL Lacs). Their spectra 



acterized by a two-bump structure (Fossati et al. 



ener gy distributions (SEDs: in vF v ) are char- 
1993). Since the discovery of the first extra- 



galactic TeV-photon emitter, Mkn421 (IPunch et al.lll992h . very-high-energy (VHE: E > 80 GeV) 
emission has been confirmed in more than 20 BL Lac objects. The SEDs of many of those BL 
Lacs show the peaks of the lower energy bump at UV to X-ray energies . These objects be- 
long to the sub-class known as "High-frequency peaked BL Lacs (HBLs)" (|Padovani & Giommi 
19951). Their non-thermal emission in this lower energy bump is commonly ascribed to syn- 
chrotron rad iation from relativist ic electrons, accelerated in the jet moving with relativistic bulk 
speed (e.g.. iGhisellini et al.N 1998D . The two-bump structure in SEDs of HBLs has been well ex- 
plained by synchrotron self-Compton (SSC) models, where the target photons of inverse-Compton 
scattering for the hig her energy bump are the synch rotron photons produced by the same elec- 
tron population (e.g., Costamante & Ghisellini 20021) . In this model the high energy end of the 
electron spectrum is responsible for both X-ray and VHE 7-ray emission. The observed correla- 
tions of the X-ray and VHE 7-ray fluxes during large flares of VHE 7-ray emitting HBLs (e.g., 
Takahashi et al.ll 19961 : LMaraschi et al.ll 1999c iKrawczynski et aU200lb provide strong experimental 
evidence for the SSC mechanism for HBLs . The target photons could also be produced in th e 
accretion disk (|Dermer & Schlickeisenl 1993b or in the broad-line region (e.g JSikora et al.lll994|) . 
Alternatively, the high-energy emission c an be also due to pions produced by accelerated pro- 
tons and ions and su bsequent pion dec ay (IMannheimlll993b or direct synchrotron emission from 
high-energy protons (|Aharonian 112000b . 



The results of applying emission models to the data can provide information on physical 
parameters of the jet, such as the co-moving magnetic field, the population of the accelerated elec- 
trons, the Doppler boosting factor and the s ize of the emitting region. HBLs often show strong flux 



variabi lity on time scales of less than 1 hr (|Gaidos et al.ll 1996c lAharonian et al.ll2007c lAlbert et al 



2007b|) . Hence, simultaneous multiwavelength observations over a wide energy range, covering 
in particular X-ray and VHE 7-ray bands, are essential to study the physics of these high-energy 
radiation emitters. 

Until a few years ago, simultaneous multiwavelength observations were only possible during 
flaring states due to the low sensitivity of the participating 7-ray telescopes. In the VHE 7-ray 
band, new generation of Imaging Atmospheric Cherenkov Telescopes (IACTs) such as MAGIC, 
H.E.S.S., and VERITAS, can access the energy range from below 100 GeV up to several TeV 
These instruments also allow us to detect GeV-TeV 7-ray signals within short observation times 
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(several hours) even in quiescent source states. A comparison of emission model parameters for 
several different source states may allow us to reveal the origin of the jet activity. 



The Suzaku X-ray satellite (see section |2T2l ) features the most sensitive instruments among cur- 
rent X-ray detectors for time-resolved coverage of a wide X-ray energy band, from the soft to the 
hard X-ray energies, well beyond lOkeV. Suzaku already performed observations of sever al known 



TeV HBLs and successfully obtained time-reso lved spectra up to hard X-ray energies (|Sato et al 



20081 : iTagliaferri et al.ll2008l : lReimer et al.ll2008|) . The capability to perform successful multiwave- 
length o bservations for TeV-HB Ls with MAGIC and Su zaku even in quiescent states has been 
shown in ITagliaferri et all d2008h and lReimer et all d2008l) . 



Mkn 501 (redshift z = 0.034) is categorized as an HBL and was the second established TeV 
blazar ( Ouinn et al. 1996|) . In 1997 this source went into a state of surprisingly high activity. The 
detected flux was 10 times higher than that of the Crab Nebula in the VH E 7-ray regime; the 



high energy photons were observed up to ~20 TeV (lAharonian et al. 



Bepp o SAX also observed a flaring activity in the X-ray band (ICatanese et al.l 1 19971 : IPian et al. 



200 lb. At the s ame time, 



19981 : iTavecchio et al.ll200l|) . The observed X-ray data showed an exceptionally hard spectrum 
with a synchrotron peak at (or above) ~ lOOke V This represents a s hift of at least two orders of 
magnitude with respect to previous observations. iGliozzi et al.l (|2006) organized a long-term mon- 
itoring campaign in 2004, also covering the X-ray and TeV energy bands. They confirm the pres- 
ence of a direct correlation between X-ray and VHE 7-ray emission, which appears to be stronger 
when the source is brighter. In 2005, when MAGIC observed the object between May and July, 
the source flux varied by an order of magnitude. During the two most active nights, rapid VHE 7 - 
ray flux variability with a doubling time of a few minutes was observed ( Albert et al] l2007b|) . 
Several exte nsive SEP studies based on multiwavelength observations of this object were re 



ported (e.g.,[K ataoka et al 



1999: iKrawczvnski et all 2000: S ambruna et alJ l2000: Ta vecchio et al. 



2001 ; G hisellini et al.ll2002|) . However, the multiwavelength observations that included VHE 7-ray 
and X-ray instruments were only conducted during flaring states, but no simultaneous X-ray and 
VHE 7-ray data were available for a low state of activity. 

In 2006 July, a joint multiwavelength campaign between MAGIC, Suzaku, and the Kung- 
liga Vetenskaplika Academy (KVA) optical telescopeQ was organized to observe Mkn 501. We 
succeeded in obtaining clear detections from the simultaneous observations both in the VHE 7-ray 
band by MAGIC and in the X-ray band by Suzaku. In this paper, we report the observational results 
of this campaign. The observations and data reduction for both instruments are briefly described 
in Section 2. In Section 3, we present the results of the measured light curves and spectra for each 
energy band. In Section 4, we discuss the application of a simple one-zone SSC model to the SEDs 



1 http : //tur3 . tur. iac.es/ 
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obtained in this campaign and compare those to the historical data from flaring states. Finally, we 
summarize our results in Section 5. 



2. Observations and data reduction 



2.1. VHE 7-Ray band: MAGIC 

The MAGIC telescope is an IACT with a 17m diameter dish, located at the Canary Island La 
Palma(28°.2N, 17°. 8 W, 2225 ma.s.L). 

In 2006, Mkn 501 was observed by MAGI C between July an d September. The observations 
were performed in the so-called wobble mode (|Daum et al.lll997|) . where the object is observed 
with an 0.4° offset from the camera center. With this observation mode, ON- and OFF-data sam- 
ples can be extracted from the same observation run; in our case, we used three OFF regions to 
estimate the background. As a part of the multiwavelength campaign, intensive observations were 
conducted during three nights on July 18th, 19th, and 20th with a total observation time of 10.5 nr. 
In eight additional nights, MAGIC pointed at this source for only a few tens of minutes each time 
as a part of an extended monitoring program. In total, 4.2 hr were spent for these additional mon- 
itoring observations. After rejecting the data with anomalous trigger rate due to bad observation 
conditions, the remaining good quality data were further analyzed. Data taken under large zenith 
angles (> 35°) were also excluded to maintain the low energy threshold. 

A shower image cleaning wa s applied based on the charge amplitude and time information 
in each pixel ( Albert et al. 2008a|) . Every cleaned event was parameterized using the so-called 
Hillas parameterization ( Hillasl 1985 ). These parameters were used for 7/hadron se paration and 



energy estimation of 7-ray events by means of the "Random Forest (RF)" method (Albert et al 



2008b|) . In the RF method, Monte-Carlo simulated 7-ray samples (M ajumdar et al. 



2005) with 



the same zenith angle range as the data were used as 7-ray training sample while real data were 
used as hadron-event background sample. The 7-ray signal was extracted on the basis of the DISP 
method (iDomingo-Santamarfa et al.ll2005|) : a cut on the 9 2 parameter (the squared angular dis- 
tance between the nominal s ource position and the reconstructed 7-ray direction) was applied to 
determine the 7-ray signal d Albert et al.l 12007 al) . Final spectra were derived using an unfolding 
technique (|Albert et al.ll2007c|) . More detailed infor mation on the st andard analysis steps and per- 
formance of the MAGIC telescope are given in Albert et al.l (|2008al) . As quoted in that paper, we 
estimate a systematic energy scale error of 16%, a systematic error of 11% on the flux normaliza- 
tion (without the energy scale error), and a systematic slope error of ±0.2. 

For the multiwavelength campaign, based on 9.1 hrs of good quality data an excess of 1513 
events over 261 12 normalized background events yielding a significance of 8.0 a was obtained for 
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the following analysis. We note that tighter cuts that only selected data of a shower image size > 
350 photoelectrons (corresponding to a 7-ray energy peak of about 250 GeV) with a 9 2 < 0.03 deg 2 
resulted in an increased 13.4 a significance. 



2.2. X-Ray band: Suzaku 



The joint Japanese-US satellite Suzaku (|Mitsuda et al.ll2007|) . launched successfully into orbit 
on 200 5 July 10, covers a wide energy range of 0.2-600 keV. Suzaku carries four sets of X-ray tele- 
scopes (ISerlemitsos et al.ll2007 ) each with a focal plane X-ray CCD camera (X-ray Imaging Spec- 



trometer (XIS); Koyama et al 



2007), covering an energy range of 0.3-12 keV. Three of the XIS 



(XIS 0,2,3) detectors have front-illuminated (FI) CCDs, while the XIS- 1 utilizes a back-illuminated 
(BI) CCD. The merit of the BI CCD is its improved sensitivity in the soft X-ray en ergy band below 



1 keV . Suzaku also features a non-imaging collimated Hard X-ray Detector (HXD; lTakahashi et al. 



20071) . consisting of PIN silicon diodes for the lower energy band (10-70 keV) and GSO scin- 
tillators for the higher energy band (40-600 keV). Suzaku observations can be conducted using 
two default pointing positions, the XIS nominal position and the HXD nominal position. In this 
observation for the multiwavelength campaign, we used the HXD nominal position to maximize 
the effective area of the HXD. In the following analysis, the HXD/GSO data are not used because 
there is no significant detection above the 3a level. 

X-ray observations by Suzaku were carried out between 2006 July 18, 18:33:00 UTC and 
2006 July 19, 17:27:00 UTC (sequence number 071727010). All XIS sensors were operated with 
1/4 window option in order to reduce possible pile-up effects. In total 35 ks of good time intervals 
(GTIs) is obtained for each XIS and HXD detector after screening criteria as described in the 
following. 

The XIS data used in this paper were reduced via the Suzaku software version 2.0. The screen- 
ing was based on the following criteria: (1) ASCA-grade 0,2,3,4 and 6 events were accumulated, 
and the CLEANSIS script was used to remove hot or flickering pixels. (2) The time interval after 
the passage of the South Atlantic Anomaly is greater than 500 s. (3) Data were selected to be 5° 
in elevation above the rim of the Earth (ELV) (20° above the day-Earth rim). The XIS events were 
extracted from a circular region with a radius of 4.2' centered on the source peak, whereas the 
background was accumulated in an annulus with inner and outer radii 5.6' and 11.1', respectively. 
The response (RMF) and auxiliary files (ARF) are pro duced using the anal ysis tools XISRMFGEN 
and XISSIMARFGEN developed by the Suzaku team (|Ishisaki et al.ll2007|) . which are included in 
the software package HEAsoft version 6.5. 



The HXD/PIN data (ver.2.0) were processed with basically the same screening criteria as 
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those for the XIS, except for ELV> 5° through night and day and a cutoff rigidity >8 GeV/c. 
The HXD /PIN instrumental background spectra were provided by the HXD team for each ob- 
servation (|Fukazawa et al.ll2006l : iKokubun et al.ll2007|) . The HXD/PIN data also include another 
background compon ent, the so-called, cosmic X-ray background (CXB). In our analysis, we use 
the CXB spectrum (IGruber et all 19991) as 



— = 9.0x 10- 
dE 



—) 

3keVy 



-0.29 



exp f- 



-E \ 



erg 



(1) 



v 40keV ) cm 2 s srkeV 

The observed spectrum was derived assuming the PIN-detector response is isotropic for diffuse 
emission. Both the source and background spectra were made with identical GTIs, and the expo- 
sure was corrected for detector deadtime of 6.0%. We use the response file version ae_hxd_pinhxdnom2_2 

Spectral analysis in the X-ray band was performed using XSPEC version 11.3.2. Each XIS 
spectrum is binned such that each bin contains at least 40 counts. After the binning, we ignored 
bins with energies below 0.6 keV and above 10 keV. We also excluded bins between 1.7 and 
1.9 keV because there exist large syste matic uncertainties in the response matrices. Based on a 
contemporaneous fit to the Crab spectra. ISerlemitsos et al.1 (|2007l) reported that the normalizations 
among the CCD sensors are slightly different (by a few %). A re lative normalization of HXD/PIN 
detector to the CCD cameras also needs to be taken into account ( Kokubun et al.ll2007 ). Therefore, 
the XIS-2, XIS-3 and HXD/PIN spectra were scaled by a constant factor with respect to the XIS-0 
spectrum. 



3. Results 
3.1. Light curves 

Figure \T\ shows the diurnal light curves for Mkn 501 in 2006 with the VHE 7-ray emission 
above 200 GeV as measured by MAGIC, public X-ray data taken by RXTEI ASM§ and optical R- 
band data provided by the Tuorla Observatory Blazar Monitoring Program^ using KVA. In 2006, 
the source generally showed a low state of activity in the VHE 7-r ay band, in contrast to 2005 



when the flux varied by an order of magnitude up to 3.5 Crab units (Al bert et al.H2007b|) . Such a 
strong flare activity is not found neither in the X-ray nor the optical bands between 2006 July and 
September. 



2 http://xte.mit.edu/ASM_lc.html 

3 more information at http://users.utu.fi/kani/lm/ 
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The light curves during the multiwavelength campaign in different energy bands taken by 
MAGIC, Suzaku and KVA are shown in Figure [2] The binning of the VHE 7-ray data is 1 hr. The 
average integrated flux above 200 GeV is (4.6±0.4) x 10" 11 cm" 2 s" 1 (x 2 /d.o.f. = 10.1/ 10), which 



corre sponds to about 23% of the Crab Nebula flux as measured by the MAGIC telescope ([Albert et al 



2007b|) . No significant variability is found. However, due to the low source flux level, we could 



only have seen variability if the flux were to increase by a factor of 2-3. 

X-ray count rates of Mkn 501 recorded with the 4 XIS detectors (0.5-10 keV) and the HXD/PIN 
detector (12-60 keV) are plotted in the middle panel of Figure |2] Each point represents the rate in 
a time interval of 1440 s for the XIS and 5760 s interval for the HXD/PIN. The X-ray count rate 
in the XIS shows clear variability (x 2 /d.o.f. = 2 x 10 4 /41 for a constant flux fit). It gradually rises 
during the observation and an overall increase of about 50% can be seen. The HXD/PIN count rate 
also seems to follow that increasing trend, but a fit with a constant value yields x 2 /d.o.f. = 13.3/14. 
The optical i?-band flux is also consistent with a constant value with an average of 16.6±0.1 mJy. 



3.2. Spectra 



Figure |3] shows the unfolded differential spectrum in the VHE 7-ray band, averaged over the 
three days during the campaign. It is well described by a simple power law from 80 GeV to 2 TeV 
(with x 2 /d.o.f. = 1.85/5): 



dF 
dE 



;(1.14±0.10) x 10 



-10 



0.3 TeV 



-2.79±0.12 



photons 
TeV cm 2 s 



(2) 



An independent analysis gives results in very good agreement with those quoted numbers. The 
flux level and the photon index of this spectrum are similar to those in the lowest state among the 
MAGIC results in 2005 which were presented on a night-by-night basis in Albert et al.1 (|2007b|) . 



Figure |4] shows the Mkn 501 7-ray spectrum obtaine d during this multiw avelength campaign 
together with four spectra measured by MAGIC in 200 5 ([Albert et alj|2007bh an d a spectrum by 



C AT Cherenkq v imagin g Telescope on 1997 April 16 (|Djannati-Atai et al 



1999). As described 



in lAlbert et al.1 (|2007bh . the 2005 MAGIC data were subdivided into three groups, i.e., low-, 
medium- and high-states depending on the diurnal integral flux level. The other MAGIC spec- 
trum corresponds to data taken during a strong flare on 2005 June 30. The CAT data were taken 
in 1997 during a previous multiwavelength campaign with the BeppoSAX X-ray satellite when 
the source was in a flaring state. Historically, spectra of Mkn 501 in the VHE 7-ray band have 
shown strong variability, and different features of the variability can be seen depending on the 
energy bands as shown in Figure |4] The difference in flux at ~ 1 TeV reaches almost two orders 
of magnitude, while a difference of only a factor of 2-3 in the flux can be seen around 100 GeV. 
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The spectrum gets harder as the source flux increases. The energy of the spectral peak, which we 
consider to be the inverse-Compton peak, can be seen above 100 GeV in the highest flux states (on 
1997 April 16 by CAT and on 2005 June 30 by MAGIC), while the peak seems to be located well 
below 100 GeV in the 2006 observations. 

In the 0.6-40 keV X-ray band, we perform a joint fitting with three FI-XISs (XIS-0, 2 and 
3: 0.6-10 keV) and HXD/PIN (12 -40 keV) to deriv e a spectrum. A fixed value of the galactic 



column density at 1.73 x 10 20 cm 2 (IStark et al.lll992|) is used for the galactic absorp tion. The nor- 



maliz ation of the HXD/PIN with respect to XIS-0 is fixed to be 1.15 as reported in 



Kokubun et al. 



J2007b . A simple power-law fit yields a photon index of T = 2.292 ± 0.002, but this model gives 



an unacceptable fit (a reduced Xy=l-107 with 3893 d.o.f: probability= 2.4 x 10" 4 %). Fitting with 
a broken power law, we obtain a significantly improved fit (a reduced xl = 1-050 with 3891 d.o.f: 
probability= 1.5%) with the best-fit photon indices of T x = 2.257 ± 0.004 and T 2 = 2.420 ± 0.012 
below and above the break energy Et, r k = 3-24^;}| keV, respectively. The constant factors for XIS-2 
and XIS-3 are 1 . 033 ± 0.003 and . 98 1 ± 0.003, respectively, which are within the values previ- 



ously reported in iSerlemitsos et al.l (|2007D . A double-broken power-law model does not improve 
the fit at all. Figure [5] shows background- subtracted folded count spectra of the three FI XISs and 
HXD/PIN with residuals for the broken power law, using the fixed Galactic column density. The 
best-fit parameters and associated errors are summarized in Table [T] These derived fit parameters 
are in good agreement with the results from the XIS data (0.6 -10 keV) fit alone. It suggests that 
there is no significant change (e.g., high-energy cut-off) in the spectrum between E br u and 40 keV. 
The photon index below the break (T\) clearly shows a softer value than 2. This indicates that the 
synchrotron peak is located at energies below the XIS range, i.e., lower than 0.6 keV. 

We also attempted to fit time-resolved broken power-law spectra. The temporal behavior of 
three parameters, Y\, T 2 , and the break energy, as well as the model flux between 2 and 10 keV, 
are presented in Figure [6] with a 5760 s interval. The normalization factors for the XISs were 
the same as stated above (see Tabled])- The model flux increases, similar to the XIS count rate. 
Figure [7] shows a scatter plot between the model flux between 2 and 10 keV and the photon index 
IY A spectral hardening trend by ~0. 15 can be seen in T\ as the flux increases with a correlation 
coefficient of r = -(0.75^2q). On the other hand, the other parameters do not show strong, flux- 
dependent variations. 

Although the X-ray data show some variability (increase by ~ 50 %), given the low flux 
level, it is still much less than the a factor of ~ 2 that would be needed to clearly detect the 
variability in the VHE 7-ray range. Thus, in the following section we discuss the broad band 
SED of Mkn 501 using average spectra in both VHE 7-ray and X-ray bands taken during this 
multiwavelength campaign. 
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4. Discussion 



Figure [8] shows the overall SEDs of Mkn 501 with data obtained during this multiwavelength 
campaign as well as some historical data. The flux in the VHE 7-ray band is cor rected for ab 
sorpti on by the extra-galactic background light (EBL) using the "low-IR" mo del of 



Kneiske et al. 



2007) 



(120041) . In our optical data, the host galaxy contribution (12.0 ±0.3) (mJy) (INilsson et al 
has already been subtracted. 

Assuming a uniform injection of the electro ns throughout a homogeneou s emission region, 
we applied a one-zone SSC model, developed by iTavecchio et al.1 (|1998L 1200 lr) . to our campaign 
data for estimating physical parameters of the emitting region. A spherical shape (blob) of radius R 
is adopted for the emission region, filled with a tangled magnetic field of intensity B. The electron 
distribution is described by a smoothed broken power-law energy distribution with slopes n\ from 
7min to the break energy 7 b and n 2 up to a limit of 7 max and with a normalization factor K. The 
relativistic effect is taken into account by the Doppler beaming factor 8. 



Th e radius, R, is selected to be 1 .03 x 10 15 cm, corresponding to the value reported in lAlbert et al 
(I2007bl) for the SEDs observed in 2005. 7 max is set to be 10 7 since no cut-off in the high energy 
ends of both the X-ray and the VHE 7-ray spectra is detected. 7 m ; n is fixed at 1 as a nominal value 
because this parameter does not affect the emission in the energy bands of our data. Since the data 
do not clearly indicate the positions of the synchrotron and inverse-Compton peaks, we cannot 
fully constrain the SSC parameters by precise fits to the SEDs. In addition, variability timescales 
cannot be determined for the data because of the quiescent state of the source. Therefore, our aim 
is to reproduce the observed spectral behavior and correlations within a simple unifying picture 
rather than entering in details of the jet structure. 

First, we applied the SSC model for the low state SED which is obtained during our multi- 
wavelength campaign in 2006. The one-zone SSC model can reproduce the measured X-ray and 
VHE 7-ray spectra in this low state of activity of the source as shown in Figure[8] However, it is ap- 
parent that the model in Figure [8] underestimates the flux in a low energy range between radio and 
optical. Usually, homoge neous models cannot be used to explain the low frequen cy radio emission 



(see e .g., IPian et al.lll998|) due to efficient self absorption. In pr evious studies. iKatarzvnski et al. 



(1200 1|) used an inhomogeneous conical jet model proposed by iGhisellini et al. 09851) to explain 



radiation from the low radio frequency up to the ultraviolet of Mkn 501. Here, along the same 
lines, we consider the energy range between X-ray and TeV for our one-zone SSC model. 

On the basis of model parameters for this low state, we attempted to reproduce the SED ob- 
tained during the flare on 2005 June 30, using the same SSC model. There are no simultaneous 
X-ray data other than RXTE/ASM available at that time. The measured flux by the RXTEI ASM 
(pink triangle in Figure [8]) shows a compatible level in the X-ray spectrum to those taken by Bep- 
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poSAX (green dots in Figure [8]) on 1997 April 16. In addition, the VHE 7-ray spectrum taken by 
MAGIC on 2005 June 30, was almost equivalent to the spectrum measured by CAT on 1997 April 
16, as shown in Figure HI Therefore, we used this BeppoSAX spectrum as a guide for the X-ray 
spectrum during the VHE 7-ray flare on 2005 June 30. With this, we can reproduce the SED in 
this high state just by changing 7 b , the Lorentz factor of the electrons at the break energy in the 
electron spectrum. The SSC models for the low and the high states of Mkn 501 are represented by 
dashed lines in Figure [8j The derived parameters for these SSC models are listed in Tabled 

In Table [3] we compare our results to some of the previous SED studies based on SSC models 
for Mkn 501. All of those were derived from applying one-zone SSC models to actual observa- 
tional data. Not all studies used simultaneous X-ray and VHE 7-ray data. In fact, most simul- 
taneous X-ray and VHE 7-ray data were taken during the huge outbursts in 1997. Nevertheless, 
the SSC model parameters of 8 and B from our multiwavelength campaign in 200 6 indicate values 
simil ar to those of previous works for different flux states, apar t from models o f 



Sambruna et al. 



hood) (lower B) for the medium flux state in 19980, and of iKonopelko et al.l feooj) (higher S 



an d lower B) for the high flux state in 1997. Those parameters are also consistent with values 



111 



Bednarek & Protherod (|1999r) . who constrain the parameters of the emission region based on 



the variability time scale during the 1997 April 15-16 flaring activity. Note again that only our 
work used simultaneous X-ray and VHE 7-ray data taken in a low flux state. 



Tavecchio et al.l (|2001|) attempted to model different emission states of Mkn 501 in 1997 and 



1999 by mainly changing the break energy of electrons, slightly modi fying its spectral slopes 
and number density, and by keeping other parameter unchanged. Also IPian et al.l (|1998i) could 
reproduce different flux states by changing only the electron distribution, keeping the same values 
for other parameters. In these frameworks, the electron spectrum is the key component representing 
the different activity states of Mkn 501 ; especially, 7b rea k can play a major role there. 

From the ph ysical point of view, in the contex t of the widely discussed diffusive shock accel- 
eration models (Kirk et al. Il998l : iHenri et alll999|) the variations of 7b rea k (and the other parameters 
specifying the particle acceleration) could be explained by changes of the parameters determining 
the efficiency of the acceleration mechanisms (such as the parameters characterizing the turbu- 
lence). A deeper discussion of this point is clearly beyond the scope of this paper. 



4 however, they did not take into account the absorption by EBL in the VHE 7-ray data. 
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5. Summary 

We present first results of a multiwavelength campaign for the TeV blazar Mkn 501 during its 
low state of activity observed by MAGIC and Suzaku in 2006 July. 

• VHE 7-ray signals were clearly detected at a 13.4 a level from 9.1 hrs of data taken during 
the campaign. The average integrated flux above 200 GeV corresponds to about 23% of 
the Crab Nebula flux without significant variability detected in the data. The spectrum in 
the VHE 7-ray band is well described by a simple power law from 80 GeV to 2 TeV with a 
photon index of 2.8 ± 0. 1 . The flux level and the photon index of this measured spectrum are 
compatible with those found in the lowest state among the MAGIC Mkn 501 observations 
in 2005. 

• The X-ray spectral shape derived from the Suzaku data from 0.6 keV up to 40 keV is well 
described by a broken power-law without a need for a high energy cut-off. The derived 
photon index suggests that the synchrotron peak is located below 0.6 keV. The flux level in 
X-ray showed a low state of activity, in similarity to the VHE 7-ray flux. As the flux shows a 
small increase (~50 %), the spectral index below the break energy shows a hardening trend 
by ~0.15. 

• The overall SED in the low state during our multiwavelength campaign can be well repre- 
sented by a homogeneous one-zone SSC model. Based on the SSC parameters for this low 
state, the high state SED in 2005 can be reproduced by changing only the Lorentz factor 
of the electrons corresponding to the break energy of the electron spectrum. This suggests 
that the variation of the injected electron population in the jet is responsible for the observed 
variation of the SED. In particular, the electron energy at the spectral break could be a key 
parameter to represent the different activity states of Mkn 501 . 
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Fig. 1. — Diurnal multiwavelength light curves during the MAGIC observations of Mkn 501 in 
2006 July-September. The vertical band represents the window of the Suzaku pointing. [Top]: 
VHE 7-ray flux above 200 GeV as measured by MAGIC. The horizontal dotted line represents the 
half flux level of the Crab Nebula in this energy range. [Middle]: averaged daily X-ray count from 
RXTE/ASM. [Bottom]: optical i?-band flux by KVA. 
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Fig. 2. — Light curves in different energy bands during this campaign. Each dotted horizontal 
line represents the average flux for each measurement. [Top]: VHE 7-ray flux measured by the 
MAGIC telescope. [Middle]: X-ray count rates measured by Suzaku with the four XIS detectors 
(filled circle) and the HXD/PIN detector (open square). [Bottom]: optical R -band flux measured 
by KVA. 
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Fig. 3. — Measured differential energy spectrum at VHE 7-rays for Mkn 501 with the MAGIC 
telescope. The blue line corresponds to a simple power law fit . The fit parameters are listed in the 
figure. For comparison, the measured MAGIC Crab spectrum ([Albert et al]|2008ah is shown as red 
dashed line. Vertical bars denote the 1 a statistical error. Horizontal bars represent the size of the 
energy bins. 
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VHE y-ray spectra of Mkn501 
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Fig. 4. — Meas ured VHE 7-ray spectra o f Mkn 501 in different activi ty states. The CAT data 
were taken from bjannati-Atai et all dl999h . the 2005 MAGIC data from I Albert et all d2007bh . the 
2006 MAGIC data from this work. Vertical bars denote the 1 a statistical error. Horizontal bars 
represent the size of the energy bins. 
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Fig. 5.— Suzaku (XIS-0 [black], XIS-2[red], XIS-3[green] and HXD/PIN [blue]) averaged spec- 
trum of Mkn 501. The model plotted with the data is a broken power law obtained by a joint fitting 
to these three XISs and HXD/PIN data. The parameters are shown in Table CD The lower panels 
show the residuals for this broken power-law model. Vertical bars denote the 1 a statistical error. 
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Fig. 6. — Temporal behavior of the fitting parameters to the XIS data with a broken power-law 
model. Each point represents a 5760s interval. Photon indices below (T\) and above the break 
energy (fit,*) an d model flux between 2 and 10 keV are described. A horizontal dotted line in each 
panel represents an average value of each parameter. Note that because of a correlation between 
T 2 and E hrk , the errors of T 2 can also vary depending on the uncertainties of E hrk . 
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Fig. 7. — Scatter plot of X-ray flux between 2 and 10 keV and photon index below the break energy 
(Ti), which are taken from Figured 
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Fig. 8. — Overall SED of Mkn 501 as measured in July 2006 and historical data. Red points rep- 
resent energy fluxes from this campaign obtained by KVA (optical), Suzaku (X-ray) and MAGIC 
(VHE 7-ray). VHE 7-ray fluxes are corrected by the "low-IR" EBL model of lKneiske et al.l(|2004l) . 
Corresponding measured VHE 7-ray fluxes are also plotted by blue points. Green points de- 
scribe some historical X-ray and VHE 7-ray fluxes. The X-ray spectra were obtained by Bep- 
poSA X on 1997 April 16 (th e highest), 1997 April 29 (medium) and 1999 June (the lowest) (taken 
from iTavecchio etaD 120011)1. The X - ray flux measured by RXTE/ ASM on 2005 Jun 30 is also 



Albert et al.1 (12 007b)). The historical VHE 7-ray spectra were obtained 



Albert et al.1 (|2007b|) . They are also corrected for EBL absorption 



shown as pink point (see 
by MAGIC in 2005 taken from 
using the same Kneiske model. Grey points a nd a bow-tie are hist orical data taken from NASA Ex- 
tragalactic database (radio-optical) and from lKataoka et al.l (11999) (7-ray data observed by EGRET 
in 1996), respectively. The l ong-dashed and short-dashed lines describe the SSC model based on 



Tavecchio et al.l (1 1 99 8l 120011) for this campaign data and the "high" state, respectively. The model 



parameters can be seen in Tabled 



Table 1 . Best- Fit parameters for Suzaku data. 



Model 
(1) 


(2) 


£brk 

(3) 


r 2 

(4) 


const. const. 
[XIS2] (5) [XIS3] (6) 


const. 
[HXD/PIN] (7) 


-Fi-lOkeV 

(8) 


X 2 /d.o.f. 
(9) 


Suzaku XIS + HXD/PIN 


power law 
broken power law 


2.292 ±0.002 
2.257 ±0.004 


i 24+0. 13 
J - -0.12 


2.420 ±0.012 


1.033 ±0.003 0.981 ±0.003 
1.033 ±0.003 0.981 ±0.003 


1.15 (fixed) 
1.15 (fixed) 


5.370 ±0.012 
5 3H+0.014 

-0.019 


1.107/3893 (2.4 X 10^% 
1.050/3891 (1.5%) 


Galactic column density: Nh = 1.73 x 10 20 [cm 2 ] 


(fixed) 











Note. — Col. (1): model used to fit the data. Col. (2): photon index for the power law model, or low-energy photon index for the broken power law model. Col. (3): 
break energy [keV] for the broken power law model. Col. (4): high-energy photon index for the broken power law model. Col. (5,6,7): constant factors with respect to 
XIS1, for XIS2,XIS3 and HXD/PIN, respectively. Col. (8): flux in the 2-10 keV band, in units of 10-" [erg cm" 2 s -1 ]. Col. (9): reduced x 2 /degrees of freedom and 
corresponding probability. 
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Table 2. SSC model parameters of Mkn 501 . 



data 


R 


" Tmin 


7br 


Tmax 


B 


K 


m 


n 2 




cm 








Gauss 


particle/cm 3 






2006 (low) 


1.03 x 10 15 


20 1 


6x 10 4 


1 x 10 7 


0.313 


8.5 x 10 4 


2 


3.9 


2005 (flare) 


1.03 x 10 15 


20 1 


1.3 x 10 6 


1 x 10 7 


0.313 


8.5 x 10 4 


2 


3.9 



Table 3. Comparison of the SSC model parameters, 8, B and R to previous studies for Mkn 501. 



5 


B[G] 


R [cm] 


Flux State 1 


reference 


15 


0.8 


5 x 10 15 


H (1997), M(1997) and L 


Pian et al. (19981 


15 


0.2 


4.5 x 10 15 


L 


Kataoka et al. (19991 


25 


0.1 


4x 10 16 


H(1998) 


Sambruna et al. (2000) 


25 


0.03 


2x 10 15 


M (1998) 


Sambruna et al. (2000") 


14 


0.15 


2.9 x 10 15 


H (1997) 


Katarzvnski et al. (2001) 


14 


0.15 


4.2 x 10 15 


M (1997) 


Katarzvnski et al. (200 1) 


10 


0.32 


1.9 x 10 15 


H (1997), M (1997) and L 


Tavecchio et al. (2001) 


11 


0.2 


10 x 10 15 


L 


Kino et al. (2002) 


50 


0.04 


3.5 x 10 15 


H (1997) 


KonoDelko et al. (2003) 


20 


0.313 


1.03 x 10 15 


H (2005), L (2006) 


this work 



'H: high state, M: medium state, L: low state (among historical data in X-ray. BeppoSAX (green) 
data in Figure [8] represent those states). Numbers represent the year when the corresponding data 
were taken. Previous works for "L" use data taken in different years in different energy bands. 



